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February 5, 2008 

To Yukiko, who never questioned the beauty of nanomagnetism 

Abstract. We calculate the response of a molecular magnet subject to a time-varying magnetic field 
and coupled to a heat bath. We propose that observations of calculated oscillations transverse to the 
field direction may be an effective way of demonstrating quantum tunneling and thus probing the de- 
tails of level repulsion. The effective model of a triangle of Heisenberg spins and weak anisotropics, 
as has been used to model the molecular magnets {V15} and {Cu^}, is used to illustrate this. 

PACS. 75.50.Xx Molecular magnets - 75.60.Ej Magnetization curves, hysteresis, Barkhausen and 
related effects - 74.25.Nf Response to electromagnetic fields (nuclear magnetic resonance, surface 
impedance, etc.) 



1 Introduction 

Studies of the quantum behavior of open systems are al- 
most as old as quantum mechanics itself. The related phe- 
nomena of relaxation and decoherence have passed from 
conceptual to pratical questions with the improvements of 
experimental technique during the last decades. One can 
now build and manipulate physical objects at the border- 
line between quantum and classical worlds. In the future, 
materials structured on such a scale may be of use in ap- 
plications such as very tiny components of magnetic mem- 
ory and quantum computers. This would require detailed 
understanding of switching and dissipation between dif- 
ferent quantum states. For such memories we can distin- 
guish "dynamical" aspects, related to writing (or reading) 
of different states, and "relaxational" aspects, concerning 
the storage itself. If such memories were to be part of fu- 
ture quantum computers, the phase of the wavefunctions 
is essential and the loss of phase central to the questions 
of the feasibility of building such a device. It is therefore 
crucial to understand and control the coupling of the small 
magnetic object to its environment. 

There is a large literature on the subject of two level 
systems in contact with an environment Tj but in the 
magnetic systems we shall consider, there are aspects not 
present in, for example, "qubits" based on Josephson cou- 
plings. In particular there is the complexity of the mag- 
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netic spectrum: as we shall see, it is essential to take into 
account more than two levels to understand the dynam- 
ics. Furthermore the nature of the environment is unclear 
and new experiments should be designed that may shed 
light on the issue. This paper will address a multi-level 
problem, calculating the response of a relatively simple 
magnetic system to a fast-varying magnetic field, includ- 
ing dissipation to a a simple model environment. In the 
context of quantum computing this corresponds to the 
"writing" aspect on a potential quantum memory storage 
device. This will be illustrated using the example of a tri- 
angle of magnetic ions, as a model of low spin molecular 
magnets such as {V15} or {Cus}. We choose the triangu- 
lar geometry because of recent experiments, to be noted 
later, and because the spectrum is sufficiently rich that it 
displays at least some of the features not seen in simple 
two-level systems. 

Experimentally, the development of new magnetic coils 
allowing for the generation of high magnetic fields vary- 
ing on very short time scales has lead to experiments that 
can probe such interesting aspects of statistical mechanics. 
We shall argue that non- equilibrium measurements with 
pulsed magnetic fields can be more informative, if prop- 
erly analysed, than equilibrium measurements for fine de- 
tails of interactions between the spins. Such experiments 
usually involve solids, where the molecular units repeat, 
and are assumed to have sufficiently weak interactions 
that they may be considered as independent, to a first 
approximation. They are therefore analysed in terms of 
a static Hamiltonian for the individual molecule, with an 
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additional time-dependent external field. The dissipation 
is induced by an environment, which represents contribu- 
tions from coupling to phonons, dipole interactions with 
other molecules, hyperfine coupling to nuclear spins. 

2 Molecular magnets and Landau-Zener 
regime 

If the interactions between spins were perfectly isotropic, 
transitions in magnetization would occur by the cross- 
ing of eigenvalues of different spins. In real magnets an- 
isotropics mix the levels and, on some scale, crossings are 
avoided. The dynamics will therefore depend strongly on 
the fine structure of the molecular levels of the spins. If 
the level repulsion is weak, as is the case if the anisotropics 
are small, it may be possible to enter a regime where 
the temperature is sufficiently low, and the coupling to 
the environment sufficiently weak, that we can be in a 
regime of pure quantum tunneling. It is for this reason 
that there has been great interest in the Landau-Zener 
tunelling regime [5] . 

If we consider only two levels with an avoided crossing, 
the transitions between the two levels, when varying a 
control parameter such as the external magnetic field, is 
described by the Landau-Zener formula [3]: 

where the "velocity" v is defined hy v = '^^"^ ^7^^"* ' ^'^'^ 
and E2 are the eigenvalues of the two levels considered. 
V is the asymptotic probability for passing to the upper 
level given the initial condition of a state in the lower. We 
note that this formula is only for the asymptotic proba- 
bility of transfer and is for zero temperature; thus it does 
not fully determine the dynamics. It does, however, suffice 
to see when quantum tunneling can be relevant. For the 
argument of the exponential to be sufficiently small that 
the probability of staying in the lowest level is other than 
unity, the "velocity" v must usually be extremely large, 
unless the gap is very tiny: e.g. for A = IK, quantum tun- 
neling is relevant for velocities greater than v f» lO^^Ts^^. 
Thus to attain this regime either A must be extremely 
tiny or v very large. Two recent sets of experiments do, 
however, approach the tunneling regime: 

(i) In the case of high spin magnets such as {Mni2} 
[5] or {Fes} Ui the two nearly degenerate states corre- 
spond to two states of total magnetization Sz = ±10 and 
tunnelling is via a series of virtual transitions to states of 
smaller \Sz \ in the multiplet of total 5 = 10. The energies 
of these states are determined primarily by the easy axis 
anisotropy DS'^.[5] Since the splitting is tiny, the tunnel- 
ing probability can be appreciable for very modest changes 
in the external magnetic field. The experiments are, of 
course, carried out at temperatures much greater than the 
energy splitting. 

(ii) In low spin magnets such as {V15} ^ or {Cua} 
[7], the active degrees of freedom can be taken to be es- 
sentially coupled spin 1/2 moments, seated on the vertex 



of equilateral triangles. In {V15} these may be thought of 
as the three spin 1/2 V^"*" ions sandwiched between two 
planes of six ions. As de Raedt et al have shown ^8], the 
extra Vanadium ions renormalize the effective interactions 
in the plane, but for the fields of interest the out-of-plane 
spins can otherwise be considered as being locked in inert 
singlet states. The gaps are generated by Dzalyoshinsky- 
Moriya (DM) interactions and are typically 10~^ K so that 
the velocity needed for Landau-Zener tunnelling would be 
10'^T.s-\ within range of recent "single-turn" magnets [9] . 
The advantage of these compounds compared to those of 
(i) is that experiments can be conducted at temperatures 
comparable to the splittings, thus a different regime of 
non-equilibrium phenomena can be studied. Furthermore 
Choi et al. have suggested that the extra degeneracy of 
different chiral states of the triangular structures may be 
useful for applications [7]- 

A problem in the interpretation of past experiments is 
that it is not necessarily clear whether one is really in the 
tunneling regime. One approach has to be to look for the 
characteristic dependence on the sweeping velocity of the 
magnetic field and compare to the predictions of Landau- 
Zener theory. The difficulty is that features found such as 
the magnetic Foehn effect 10,, where extra plateaux in 
the magnetization occur in the non-equilibrium response, 
may be described either as phonon bottleneck effects [TT| 
[T2|IT3]. i.e. that the bath does not relax to equilibrium on 
the timescale of the experiment, or by a dynamical effect 
of the system coupled to a bath that stays in equilibrium. 
In what follows we shall look for extra features which could 
distinguish effects intrinsic to the small quantum system 
from such bottleneck effects. 



3 Model for the molecular magnets {V15} or 

{CU3} 

In order to describe the reaction of the magnet to the 
external magnetic field, we have to precise the internal 
Hamiltonian for the molecular magnets. For systems (ii), 
to which the rest of the paper is devoted, we can write the 
Hamiltonian for three spin ^ moments 

3 

i=l 

The system has 2'^ = 8 energy levels. The exchanges are 
chosen for triangular symmetry, i.e. the isotropic exchanges 
are equal and the three vectors defining DM interactions 
are determined by a single vector D and rotated around 
the three-fold symmetry axis [8] . Before discussing the dy- 
namics, we shall discuss the static properties of Tig- In 
Figure [1] we show its spectrum , i.e. with H{t) = HzBz 
constant for different values of Hz- At large fields the 
spin I components have lowest energies and cross with 
the two degenerate spin i levels of opposite chirality. To 
make the degeneracies apparent a small chiral-breaking 
term has been added. We choose this additional term to 
be a Sen-Chitra interaction Hchiiai ~ Achirai (H ■ Bz) Si ■ 
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(a) Energy spectrum as a function of the magnetic field. Both 
energies and fields are in units of J. 




Field (Units of J) 



(b) Equilibrium magnetization (in units of /ib) in the direction 
of the field for T = O.IJ. 

Fig. 1. (Color online) Spectrum of Hamiltonian TLq and equi- 
librium magnetization. 

{S2 X 5^3) [131 . The small Dzalyosliinsky-Moriya interac- 
tions give avoided level crossings. Below the spectrum 
we show the equilibrium magnetization at temperature 
fceT = O.IJ. It is clear that unless extremely low temper- 
atures are available, little information on the details of the 
level crossing can be found, other than the points where 
levels cross, i.e. the value of the isotropic exchange. 

In contrast to the relatively featureless equilibrium 
magnetization it will be seen that non- equilibrium ob- 
servables such as magnetization are much richer and we 
now turn to the question of dynamics. We remark also 
that in such triangular-based spin systems, switching be- 
tween different states involves more than just two states 
at a time; because of the chirality the crossings involve 
three (between spin | and the two chiral ^ states) or four 
(between the degenerate chiral doublet at = 0). At 



the least a generalization of the Landau-Zener form is re- 
quired. We remark that in certain cases the asymptotic 
transition probabilities for such multi-crossing may be re- 
duced, via an "Independent Crossing Approximation" |15j 
to the calculation of two-level crossings. Because of the 
importance of avoided crossings we need a general theory 
that applies for multiple crossings with the details of the 
anisotropy [16 . For direct comparison with experiment it 
is also useful to have a general form of the time-dependent 
external field generated by the coils. 

4 Environment: model and coupling to the 
magnet 

In current experiments and probably future applications, 
the magnet can not be considered as a isolated system. It 
is indeed surrounded by an environment, responsible for 
relaxation and decoherence phenomena. In order to take 
into account those effects in our calculation, we have to 
define an environment explicitly. Unfortunately the true 
environment is not yet known precisely. We therefore focus 
here on a simple environment [17j modelled as an infinite 
set of bosonic harmonic oscillators: 

Henv = duj UjI{uj) (^blb^ + ^ , (2) 

We choose the spectral density I{lu) to be a power law, 

I{uj) =cu°'0{uj), (3) 

with a — 2 (Super-Ohmic regime, as appropriate for a 
three dimensional phonon bath) , and where is the Heav- 
iside distribution. Several additional hypotheses will be 
needed in the following. First we will supposed that the 
initial state of environment is a thermal state, that is to 
say that the environment is in thermal equilibirum with a 
thermostat at a given temperature T. Furthermore the en- 
vironment is assumed to be Markovian, i.e. its correlation 
time is much shorter than the time scale of the dynam- 
ics. Taking these two hypotheses together implies a static 
environment, described at all times by the density matrix 

Pcnv (i) oc exp (-/?77e„v) Emm]. 

The coupling of the spin system to the bath is included 
with the term 

/-I-CX3 
duj I{Lu){bl+b^)X (4) 
-oo 

The operator X is some combination of spin operators Si 
and couples the spin system linearly to the boson opera- 
tors of the bath, bui bl^, with dimensionless strength g. It 
has sufhciently low symmetry that it can restore equilib- 
rium from any initial state, e.g. X = i'^i'^f + ^i'^i + 
CiSf). The coefficients ai,bi,Ci differ from spin to spin, 
in order to induce transitions between states of different 
chiralities. 
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5 Bloch-Redfield approach to the reduced 
density matrix dynamics 

The system (magnet + environment) is a closed quantum 
system described by the Hamiltonian 



'Htot — Ti-oit) + He 



coupling 



(5) 



All the physical information of this closed system can be 
determined by the density matrix ptot{t)- As the experi- 
mental set-up allows measurements on the magnet only, 
and not on the environment, we are interested in this mag- 
net as an open quantum system, described by the reduced 
density matrix p{t) — Tronv [ptot(i)]- A naive approach to 
study the magnetic dynamics would be to integrate the 
Liouville equation, coming from first principles, to obtain 
the total density matrix ptotit) at all time, and then to 
perform the trace on the environment degrees of freedom. 
However such an approach is essentially impossible to fol- 
low in practice given the high dimension of the Hilbert 
space. An approach to this problem, that is by now stan- 
dard, is to make the additional hypotheses on the envi- 
ronment presented in the previous section, to obtain a 
dynamical master equation for the reduced density ma- 
trix. This perturbative approach dates back to the 1950's 
[20 l [2T l f22] . We would like to stress that the equation does 
not come directly from first principles: it relies on addi- 
tional assumptions whose validity needs to be tested for 
each particular system. For systems with a larger number 
of levels, the equations have been put in a compact form 
[8ll24] to give the equation for the time dependence (in 
units where h = 1 )of the density matrix projected onto 
the Hilbert space p{t) of the molecular system (in our case 
the three coupled spins): 

1^ = -z[Ho(i), Pit)] - 9{[X, Rp{t)] + [X, Rp{t)]^) . (6) 
If we define the instantaneous eigenbasis {\k,t)} of Hoit) 

by 

Hoit)\k,t)^Ek{t)\k,t), (7) 
the operator R{t) is determined by its matrix elements 

(fc, t\R\m, t) = \l{ukm{t)) oJkmit))] (8) 

where uikm{t) = Ek{t) - Emit) and 

n0{u;) = {ef'--l)-\ (9) 

The external field H is taken to have arbitrary time vari- 
ation in a single direction e^, By or e^. 

We shall now present results based on integration of 
the quantum master equation ^ . In practice we integrate 
by fourth order Runge-Kutta approximation which leads 
to numerically convergent results, at least for weak cou- 
pling g and for times long enough to define several hys- 
teresis loops. The hamiltonian must be diagonalised for 
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Fig. 2. (Color online) Magnetization as a function of time. The 
applied field is H{t) = H4t)e^ 1^. H,{t) = -Hocos{nt). 
The blue lines shows the equilibrium magnetization for each 
instantaneous value of the field. The parameters used are the 
following: Ho = J, O = 2ty x IQ-^J/ft, = 0.05J, Achirai = 
J/20 and T = 0.2J. The coupling 3 = 4 x 10"^ 

each time step in order to define the matrix R{t). The 
time variation of the external field can be arbitrary but 
for our purposes we shall take it as cosinusoidal: 



H{t) = HQCOs{m)e, 



X, y, 



(10) 



and as initial conditions p{t = 0) = exp(— /37io(i = 0)), 
i.e. the spin system in equilibrium with the heat bath at 
inverse temperature (3 — -ff-rp- 



6 Numerical observations 

We shall now make calculations of the observables for a 
model case in which we take a single component of the 
DM vector along the z-direction, D = DzBz- The purpose 
here is not to fit a particular experimental system, more 
to make general statements. We first show the three com- 
ponents of magnetization as a function of time. The initial 
condition is the equilibrium density matrix for the same 
temperature and the field H{t = 0) = — J. It is seen that 
at the relatively weak couplings that can be integrated 
numerically (5 = lO^'^ to 10"**) the magnetization "lags" 
the equilibrium magnetization. After one full oscillation 
the profile (in time) of the out-of-equilibrium magneti- 
zation approximately repeats, converging slowly towards 
some steady-state form. We remark that the period for the 
sinusoidal external field is 10'' in units of J^^ . Returning 
to standard units, where Ti is not unity, this corresponds 
to a time of 5 x 10~^ seconds for J=l K. 

The information in Figure [5] can be represented in a 
form often used by experimentalists, i.e. as a hysteresis 
loop, where the magnetization is plotted a gains t the ex- 
ternal field. This we do in Figures [3(a)] and [4(a)| The dif- 
ference between the two hysteresis loops is that in Figure 
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(a) Hysteresis loops Mz{Hz). (b) Detail of the spectrum. 

Fig. 3. (Color online) Energy spectrum and dynamics of the magnetization parallel to the applied field, for H — Hz{t)ez 
parallel to Dzalyoshinsky-Moriya vector D = D^Bz- The parameters are Dz = J/20, Hq — J, fi = 2-k x 10~''j/?i, g — 10~^ 
and T = 0.2J. In the righthand panel, the details of the spectrum of the four lowest energy levels around zero field of Fig. (TJ 
showing the crossing of the two nearly degenerate doublets. 




(a) Hysteresis loops Mx{Hx)- (b) Detail of the spectrum. 

Fig. 4. (Color online) Energy spectrum and dynamics of the magnetization parallel to the applied field, for H = Hx{t)ex 
perpendicular to Dzalyoshinsky-Moriya vector D — DzBz- The parameters are identical to those of Figure O Dz = J/20, 
Ho ~ J, n = 2iT X 10"'^^, p = 4 X 10"'' and T = 0.2 J. In the right hand panel the avoided level crossing of the doublet states 
at zero field. 



3(a)[ the DM vector is parallel to th e external field direc- |4(b)| this difference is apparent. For rux {Hx) there is an 



tion and in the second, Figure 4(a) it is perpendicular to 
the DM vector. In each case the magnetization is longi- 
tudinal, i.e. parallel to the applied field. The differences 
between the two orientations {Hz) and vix {Hx) are 
that the hysteresis loop is "fatter" when the field is per- 
pendicular to the DM vector, but, more strikingly, that 
there are extra oscillations visible in nrix {Hx) that are 
absent for ruz {Hz)- What is the difference between the 
two field orientations in terms of the spectrum of excita- 
tions? If we look at an enlargement of the sp ectrum of 



Figure [T] around zero field, as shown in Figures 3(b) and 



avoided level crossing at zero field with a minimum split- 
ting whereas for rUz {Hz) there is a crossing of levels. This 
avoided level crossing makes the hysteresis loop "thinner" 
as the lowest level smoothly evolves from one spin state 
to another. So what is the significance of the oscillations? 
For a two-level system formulated as a (pseudo-) spin in 
transverse field and time-varying longitudinal field, these 
transient oscillations [24] correspond to precession around 
an instantaneous axis which, near the avoided crossing, is 
perpendicular to the applied field. In other words these 
are simply Rabi oscillations in the presence of a time- 
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Fig. 5. (Color online) Details of the oscillations in the magneti- 
sation transverse to the applied field, perpendicular to the DM 
vector. Inset: Wave packet on a longer time scale. The param- 
eters are = J/20, Hq = J , fl ^ 2-k x fO"^J/;i, 5 = 8 x lO-^* 
and T = 0.5J. The avoided crossing is centred at t=250. 



varying gap. These oscillations are not periodic, but the 
characteristic frequency range is the gap in the avoided 
level crossing. This mixing of longitudinal and transverse 
oscillation is well known in "zero-field" NMR [23] . This 
suggests that the oscillations may be more visible in the 
transverse components of magnetization, i.e. the compo- 
nents perpendicular to the applied field. In Figure [5] we 
show the transverse magnetizations, for applied field per- 
pendicular to the DM vector. It is seen that they are small 
except close to avoided crossings, where there are oscilla- 
tions rapid on the time scale of the driving field. We re- 
mark that they are visible at a temperature (0.5J) well 
above the gaps induced by anisotropics. 

This leads to the suggestion for experiment that in- 
stead of looking exclusively at the shape of the hysteresis 
loop, as has been done most in the past; it may be more 
instructive to look at the Rabi oscillations. They may be 
most visible in the transverse components, as they are 
against a flat background. We remark that these oscil- 
lating components have been seen clearly in past numer- 
ical calculations [^[TCTif^ of the longitudinal magneti- 
zation but have attracted relatively little attention. We 
recall however from the theory of two-level systems that 
the decay of longitudinal components and transverse os- 
cillations at a fixed magnetic fleld determines two char- 
acteristic time scales which give different information on 
the spectrum of the bath: in particular the transverse os- 
cillations are sensitive to the zero frequency component 
of the bath By analogy, measuring both longitudi- 
nal and transverse components one should, in principle at 
least, determine different properties of the spectrum of the 
bath. 

We remark that for quantitative studies of the oscil- 
lations we should be a little more precise about what we 
mean by "transverse" ; ideally we would take a component 



with vanishing mean. For the simple case we have shown 
here the ideal component is that which is perpendicular 
to both the applied field (parallel to e^) and parallel to 
the DM vectors, i.e. {Hx). The component perpendic- 
ular to both the field and the DM vector rUy {Hx), has 
small but non-zero values depending on the field. For non- 
coUinear DM vectors on the three spins there may be no 
absolute separation of smooth and oscillating contribu- 
tions; nevertheless the oscillations will always be stronger 
perpendicular rather than parallel to the applied field. 

We note that these oscillations occur when the rota- 
tional symmetry is broken. This occurs if, as is the case in 
Figure [51 the DM vector and the external field completely 
break rotational symmetry, so that the phase transverse is 
determined by the (time-dependent) Hamiltonian [27] . Os- 
cillations also occur, (but are not shown in the Figures), in 
much smaller amplitudes of oscillations for (Hz); that 
is even when the Hamiltonian Tio is symmetric around the 
field direction. This can happen because the coupling to 
the bath breaks the symmetry in our simulations. Whether 
this is true in a physical system or not, depends on the 
symmetries of the environment of the magnetic molecule. 
To restore symmetry in our formulation, while nonethe- 
less allowing relaxation, we could include a sum of several 
terms of the form written for Ticoupiing- 



7 Conclusions 

We have calculated the dynamic behaviour of weakly an- 
isotropic triangular antiferromagnets, as could be seen in 
the molecular magnetis {V15} or {Cua} subject to vary 
fast varying magnetic fields. Our simulations model in de- 
tail effects of tunneling including transient effects in a sit- 
uation where Landau- Zener expressions are not sufficient 
to determine dynamics at arbitrary frequency and at fi- 
nite temperature. This is especially so here are several 
levels contribute to the response because of degenerate 
chiral states on the triangle. We have argued that a use- 
ful characterization of the tunneling regime, useful would 
be the (non-periodic) oscillations in the magnetizations 
and, in particular, in the transverse components close to 
regions of level anti-crossing. This should yield more de- 
tailed information on the anisotropics causing level repul- 
sion in molecular magnets, and the coupling to the bath, 
e.g. the value of the coupling g and the nature of the bath. 
The oscillations can persist (at least at weak coupling) up 
to temperatures well above the scale of the avoided level 
crossing. To illustrate this we have calculated hysteresis 
loops for a model in which there is a single component 
to the DM vector and the field is taken either parallel or 
perpendicular to this vector. We chose the simple case of 
a single component DM vector in order to illustrate the 
phenomenon. For this case there is only a single region 
of avoided crossing at zero field. In real systems such as 
{Cus} all three components are present. There are then 
avoided crossings where, in the isotropic limit, the spin 
I cross the spin i levels H = ±| J, and oscillations will 
occur at all crossings. 
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We conclude that if one is interested in studying the 26. Cohen- Tannoudji C et al., Atom-Photon interactions, Wi- 

physics at the scale of a single molecular magnet, it might ley (1988). 

be useful to study not only the traditional hysteresis loops, 27. We thank H. Nojiri for raising this issue. 

but also the transient oscillations, visible in the transverse 
component of the magnetization dynamics, as the first cor- 
rection to the adiabatic limit. 
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